Universal Pion Freeze-out in Heavy-Ion Collisions 
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Based on an evaluation of data on pion interferometry and on particle yields at mid-rapidity, 
we propose a universal condition for thermal freeze-out of pions in heavy-ion collisions. We show 
that freeze-out occurs when the mean free path of pions A/ reaches a value of about 1 fm, which 
is much smaller than the spatial extent of the system at freeze-out. This critical mean free path is 
independent of the centrality of the collision and beam energy from AGS to RHIC. 

PACS numbers: 25.75.-q, 25.75.Gz 



A systematic study of the space-time extent and the 
dynamical behavior of the pion source in relativistic 
heavy ion collisions at thermal freeze-out can be obtained 
from analysis of pion interferometry (HBT) data. Under- 
standing these aspects is vital for interpretation of the 
data in terms of formation of the quark-gluon plasma. 
Indeed recent HBT results from RHIC and how they fit 
into the systematics have been noted as a major puzzle 
jjj . In this letter we present an investigation of the freeze- 
out conditions at beam energies from AGS to RHIC. In 
particular, the recently published CERES HBT data at 
40, 80, and 158 A GeV || provide an important link be- 
tween the existing results from AGS, SPS, and RHIC, 
thereby shedding light on the RHIC puzzle. 

Thermal freeze-out of pions and its connection to the 
mean free path has been discussed previously (see e.g. 
ii Hi H> @i □)■ The mean free path of pions at freeze-out 
is defined as 



can be replaced by the number of particles N contained 
in the freeze-out volume Vf. 

The pion freeze-out volume Vf can be accessed experi- 
mentally by pion interferometry. Mid-rapidity pion HBT 
data have been published from central collisions of lead 
and gold nuclei over a wide range of beam energies. Here, 
we focus on recent HBT results from three experiments 
which have kinematical access to the region of low trans- 
verse pair momentum k t =\\pt,i+'Pt,2\'- Experiment E895 
at the AGS §, the CERES/NA45 experiment at the 
SPS § |(J, and the STAR experiment at RHIC (TJ. All 
three experiments employ large volume Time Projection 
Chambers (TPCs), thereby applying comparable analysis 
methods with similar sources of systematic uncertainties. 

For the calculation of the freeze-out volume Vf we use 
the following expression: 



V f = (2tt)* 
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where a is the total cross section of pions with the sur- 
rounding medium and pf is the freeze-out density which 



assuming a density distribution of Gaussian shape in all 
three dimensions. The longitudinal and sideward radius 
parameters i?i on g and i? s ide are measured in the longitudi- 
nal co-moving system of the pion pair, using the cartesian 
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decomposition of the three-momentum difference vector 
<f=(<7iong,<?sidc,<7out) as proposed in Q. 

The definition of a freeze-out volume in heavy ion colli- 
sions has to be taken with some care because of the strong 
collective expansion of the system. Collective expansion 
leads to space-momentum correlations of the emitted pi- 
ons, resulting in a fcj-dependent reduction of the observed 
HBT radius parameters as compared to the 'true' geo- 
metric dimensions of the source [|l3| [l4|, [ll| . We prefer 
to use a definition containing only measured quantities 
rather than model dependent parameters that may not 
be appropriate at all beam energies. But to avoid strong 
bias of the extracted beam energy dependence due to ex- 
pansion we have selected and compared data measured 
at similar (k t ) values of about 0.16 GeV/c: The E895 
data were taken at (fe t )=0.148 GeV/c, while the STAR 
data are at (fc t )=0.170 GeV/c. From the CERES data we 
have calculated the average of the results in the fc t -bins 
at (fc 4 )=0.125 GeV/c and (fc t )=0.195 GeV/c. 

At all beam energies under investigation, the observed 
strong decrease of i?io ng as function of k t is consistent 
with a boost-invariant expansion in longitudinal direc- 
tion In this limiting case, the geometric size of the 
system in longitudinal direction is much larger than the 
measured homogeneity length i?i on g- The corresponding 
homogeneity scale in velocity space is given by the av- 
erage thermal velocity ((3th) of the pions. The thermal 
velocity ((3th) can be calculated relativistically using the 
expression: 



(ith) - 1 = - 
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with (7th)=(l - {/3th) 2 ) 2 , m t =(ml + fc t 2 )2, the ther- 
mal freeze-out temperature Tf, and Bessel functions 
Ki and K 2 - Using Eq.(||), we obtain (f3 t h)—0.7 at 
fc t =0.16 GeV/c, assuming T/=120 MeV. Turning this 
into a rapidity results in y t /j=arctanh(0.7)=0.87. We 
conclude that a source element of longitudinal size i?i ong 
corresponds to 0.87 units of rapidity at k t =0.16 GeV/c. 
Note that i?i on g is the rms of the length of homogeneity 
and so is the corresponding width in rapidity. As for the 
spatial distribution we assume a Gaussian shape also in 
rapidity. 

Also the transverse radius parameter i? s id e is reduced 
in the presence of collective transverse flow. Model de- 
pendent approximations JT7| , p"8[ |l9[ ] predict a reduc- 
tion by about 20-25% as compared to the true ge- 
ometric transverse source size at kt— 0.16 GeV/c for 
T/=120 MeV and a typical average transverse flow ve- 
locity (u t )=0.5c H |C], [H]]. We prefer not to use this 
reduction in the analysis but rather include the corre- 
sponding model dependent underestimate of the volume 
in the systematic errors below. The observed beam en- 
ergy dependence of Tf and (v t ) Q imposes an addi- 
tional systematic uncertainty. Within reasonable limits 



FIG. 1: (a): The freeze-out volume Vf as function of ^J~s. 
(b): Mid-rapidity densities of negative pions, protons, and 
antiprotons vs. yfs. 



of Tf and (vt), the relative change of i? S idc however is 
small and even partially compensated, as both Tf and 
(vt) increase with beam energy. 

The experimentally determined freeze-out volume Vf, 
calculated according to Eq. (Q) as function of ^fs, is shown 
in Fig. 0(a). We observe a steep decrease of Vf at 
AGS energies and an increase throughout the SPS en- 
ergy regime towards RHIC. The data indicate the exis- 
tence of a minimum between AGS and SPS energies. The 
origin of this non-monotonic behaviour cannot be under- 
stood in terms of previously presented freeze-out scenar- 
ios, where constant particle density in coordinate space 
or constant pion phase space density have been proposed 
as possible universal freeze-out conditions ^||, |24| . This 
is demonstrated in Fig. 0(b), where a compilation of the 
mid-rapidity densities dN/dy\ of pions and protons 
as function of ^/s is shown f^'fel M M, M, Ejfi M. 



3 



The number of pions and the sum of pions and nucleons 
are monotonically increasing with y/s, without indication 
of a minimum. 

For the determination of the critical mean free path 
A/, we need to evaluate the denominator in Eq. ([!]). In 
the presence of different particle species i we replace N-a 
by a sum: 



N ■ 



, Nia w i = N N ■ ovjv + ■ a n 



(4) 



For simplicity we neglect less abundant particle species 
such as kaons and deuterons. Following our previous ar- 
gument we integrate over the rapidity distribution with 
an rms of 2/^=0.87 and account for the different isospin 
states, leading to: 

N N = 2 ■ V2n ■ 0.87- dN p+p /dy | , and 

N n = 3-V2^-0.87-dN„-/dy I . 



(5) 

For the cross sections we use cr T7r = 13 mb, and a thermal 
and isospin averaged value a ■ 7vN =72 mb (?[. The latter as- 
sumes a thermal freeze-out temperature =120 MeV in- 
dependent of beam energy and ignoring possible changes 
due to the non-equilibrium nature of thermal freeze-out. 

In Fig. ||(a) N-a is shown as function of y/s. At AGS 
energies N-a is dominated by the contribution from nu- 
cleons while at RHIC the contribution from pions exceeds 
that from nucleons, due to the change in chemical fireball 
composition. It is the large ratio of a^N /ov^ combined 
with the steep increase in pion multiplicity with beam 
energy that produces a non-monotonic beam energy de- 
pendence of N-a, leading to a minimum in the transition 
region between AGS and SPS energies, in striking simi- 
larity with the observed behaviour of Vf in Fig. |l|(a). To 
demonstrate this quantitatively we have superimposed 
Vf and N-a in Fig. |(b) . According to Eq. (|) we conclude 
that thermal pion freeze-out occurs at constant mean free 
path of A/ « 1.0 fm. We note that A/ is much smaller 
than the system size: Thermal freeze-out is determined 
by the product of cross section times local density. A sim- 
ilar local pion freeze-out criterium was discussed in [|[ 0| 
albeit only for fireballs with comparable nucleon and pion 
numbers. This behaviour indicates a significant degree of 
opaqueness of the pion source j^], as it arises naturally 
for a collectively expanding system |32| . Considering sys- 
tematic effects on the evaluation of A/ such as the us- 
age of i?side for the volume, the Gaussian assumption for 
the spatial and rapidity distributions, the uncertainty in 
temperature (assumed as 120 MeV) and thereby in cross 
sections, and the uncertainty in evaluating the average 
thermal velocity, we find that the first two completely 
dominate, leading to 0.7 < A/ < 1.4 fm. 

Pion freeze-out at constant A / implies that the freeze- 
out density pf is constant, as long as the chemical com- 
position of the fireball does not change. This is in- 
deed demonstrated by recent results of the CERES col- 
laboration for the freeze-out volume Vf in 40, 80, and 
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FIG. 2: (a): Beam energy dependence of N-a (for explanation 
see text), (b): Comparison of N ■ a and Vf. 



158 AGeV Pb+Au collisions § shown in Fig. |(a). 
Since, at fixed energy, particle abundances scale approxi- 
mately linearly with the number of participating nucleons 
(iVpart) HJ> Hi' ^ e observation of a linear increase 
of Vf with centrality is consistent with freeze-out at con- 
stant density. 

However, as function of beam energy, relative particle 
abundances do change, and therefore the freeze-out den- 
sity pf is in general not constant. In Fig. [|(b) we show 
the pion and nucleon freeze-out densities pf^ — N^/Vf 
and pf t N=Nx/Vf as function of ^fs. The freeze-out den- 
sities change drastically with ^/s at AGS and SPS en- 
ergies, however, only little change is observed from top 
SPS energy to RHIC. This indicates that the freeze-out 
densities reach asymptotic values at high beam energy. 

In conclusion, we have derived a universal condition for 
pion freeze-out from pion interferometry data and single 
particle yields. Thermal pion freeze-out occurs at a criti- 
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FIG. 3: (a): Freeze-out volume as function of (JVpart) m 
Pb+Au collisions (from f|). (b): Pion and nucleon freeze-out 
densities pf iir and p; tN as function of y^s. At -^5=130 GeV, 
results from tt + tt + and tv~it~ interferometry are shown. 



cal mean free path A t «1 fm, independent of beam energy. 
We observe a transition from nucleon to pion dominated 
freeze-out between AGS and SPS energies, characterized 
by a minimum of the freeze-out volume Vf . The existence 
of this minimum appears as a consequence of the relative 
change of particle abundances with beam energy, under 
consideration of their different cross sections with pions. 
In this picture, the overall weak y/s dependence of HBT 
radii up to RHIC energies finds a simple interpretation. 
The surprisingly small value of A/ points to a consider- 
able opaqueness of the source || p2[ . This opaqueness 
may be at the root of the small observed values for i? ut 
at RHIC g§,g. 
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